Time independent description of rapidly oscillating potentials 
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The classical and quantum dynamics in a high frequency field are found to be described by an 
effective time independent Hamiltonian. It is calculated in a systematic expansion in the inverse of 
the frequency (u) to order co~ 4 . The work is an extension of the classical result for the Kapitza pen- 
dulum, which was calculated in the past to order uj~ 2 . The analysis makes use of an implementation 
of the method of separation of time scales and of a quantum gauge transformation in the framework 
of Floquet theory. The effective time independent Hamiltonian enables one to explore the dynamics 
in presence of rapidly oscillating fields, in the framework of theories that were developed for systems 
with time independent Hamiltonians. The results are relevant, in particular, for exploration of the 
dynamics of cold atoms. 
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The classical and quantum dynamics of a particle in a 
field that oscillates rapidly relative to the motion of the 
particle will be studied. The variation of the field in space 
is smooth but otherwise arbitrary. Such fields are applied 
experimentally to cold atoms, where a very high degree 
of control is possible. The exploration of the dynamics 
of cold atoms in strong electromagnetic fields resulted in 
many novel, interesting experimental observations 0, H 
3]. The results of this letter 0, are expected to enable 
further discoveries in this field as well as in related fields. 

For atom optics, the effect of the internal degrees of 
freedom on the center of mass motion is important. The 
force on the center of mass due to the internal degrees of 
freedom that couple to the external field is given approx- 
imately by a dipole force 0- The motion of the atoms 
can be manipulated by fields with amplitudes which vary 
spatially, resulting in a force on the center of mass of the 
atoms. In the present work, the effect of the laser on the 
center of mass motion is modeled by a time dependent 
potential. For some situations of physical interest this 
simpler model still describes the dynamics of the center 
of mass without the need to specify the dynamics of the 
internal degrees of freedom or the quantum aspects of 
the light field. In the situation that is of interest for the 
letter the frequency of this potential is large compared to 
the inverse of the characteristic time scale of the center 
of mass dynamics. An example of such systems that was 
recently realized experimentally and that motivated the 
present work is of atomic billiards HH|. The boundary of 
the billiards is generated by a laser beam that rapidly tra- 
verses a closed curve, which acts as the boundary of the 
billiard, approximated by the time average of this beam. 
The force applied by the boundary on the particles is 
approximately the mean force applied by the beam. 

The influence of a high frequency field on a classical 
particle was derived by Landau and Lifshitz |jj general- 
izing the work on the Kapitza pendulum that is a 
classical pendulum with a periodically moving point of 



suspension. The motion was separated into a slow part 
and a fast part. The leading order (in the inverse fre- 
quency) of the slow motion was calculated. This mecha- 
nism is used to trap ions in electromagnetic fields. The 
most notable example is the Paul trap 9] that can be de- 
scribed approximately by a Hamiltonian of a time depen- 
dent harmonic oscillator [ToL fill fl^ | . It is of interest to 
solve more general problems even if only approximately. 
The work of Kapitza was first extended to quantum me- 
chanical systems in a pioneering paper by Grozdanov and 
Rakovic |lj3]. They introduced a unitary gauge trans- 
formation resulting in an effective Hamiltonian that de- 
scribes the slow motion. In that paper the analysis is 
restricted to a driving potential that has a particularly 
simple time dependence. Moreover, the final results are 
restricted to forces that are uniform in space, a situa- 
tion natural in standard spectroscopy, but too restrictive 
for the interesting problems in atom optics. These re- 
strictions are avoided in the present work. Many inves- 
tigations of periodically driven quantum systems were 
performed El El El E3 El El El E3 

In the present work, a coherent theoretical treatment 
of the dynamics of classical and quantum high frequency 
driven systems is presented. Classically the motion 
can be considered approximately as one that consists of 
a rapid oscillation around a slowly varying trajectory. 
Therefore for the classical problem the motion is sepa- 
rated into a "slow" part and a "fast" part and a system- 
atic perturbation theory is developed for the motion of 
the "slow" part. This slow motion is computed to order 
and demonstrated to result from an effective Hamil- 
tonian. It is an extension of the order lu~ 2 (that is pre- 
sented in Ref. Q). Floquet theory is used to separate the 
slow and fast time scales in the corresponding quantum 
problem. An effective (time independent) Hamiltonian 
operator is defined following [13j . The eigenvalues of this 
operator are the quasienergies of the system. This ef- 
fective Hamiltonian is then computed perturbatively (to 
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order w~ 4 ). It is obtained by a gauge transformation that 
is simply related to the canonical transformation leading 
to the corresponding classical effective Hamiltonian. 
A model Hamiltonian for the motion in a periodic field 

is 

H = p 2 /2m + V Q (x) + V 1 (x,iot) (1) 
leading to Newton's equation 

mx = -Vq{x) - V((x, cot), (2) 



where V% is a 2n periodic function of tot and its average 
over a period vanishes. (We denote X = ^f, Vq = 
etc.) An instructive example of such a system is 



Vo = 0, Vi = 7 e 



-fix 1 



cos(ut). 



(3) 



The system is of particular interest since: (a) the time 
average of the potential vanishes, consequently any in- 
teresting effect is due to the rapidly oscillating potential; 
(b) when x — * oo the potential vanishes and therefore 
one expects to find scattering quasienergy states. 
We look for a solution of J5J that has the form 



,{t)=X(t)+£,(X,X,ujt) 



(4) 



where X is the slow part, X = ^ while £ is the fast 
part which is periodic (in the variable Lot) with vanishing 
average. This determines uniquely the functions X and 
£. One could use a different functional dependence of £ 
on X and its time derivatives, but such a change cannot 
affect the dependence of the solution for £ on time. Our 
method of solution is to choose £ so that J5J) leads to an 
equation for X which is explicitly time independent. An 
exact solution using is too complicated to obtain in 
general. However, at high frequencies, one can determine 
£ order by order in 1/lo, using 
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(5) 



The £j are chosen so that the equation for X does not 
depend on r = tot. One may also expand X in powers of 
1/lo as X = Xi/uj l . When one does so the equation 

of motion for X is replaced by a series of equations for 
Xi. In this series of equations, each Xj can be determined 
from the lower order terms Xj, where j < i. This is the 
standard method of separation of time scales |22j . These 
equations are equivalent, in any order, to the equation 
of motion of the (unexpanded) X which will be used in 
what follows. At a given order LO~ n of the present calcu- 
lation all contributions that are found by the method of 
separation of time scales are included, but some of the 
higher order terms are included as well. 

The slow motion is found to be controlled by the 
Hamiltonian \4i 



where g(X) = ^ \J (2)t [V{']j and P is the momen- 
tum conjugate to X. We define the integrals J^ T [/] = 



[/] 



where f(x,r)= £ n#0 fn{x)e mT and 



j times 

V if] = En^o Tnf ne " lT ■ T ne effective potential is found 
to be 



V eff (X) ee V Q (X) + V 2 (X) + Vi{X). 



where V 2 (X) 



2moj- 



■crw and 



v{' (; (2)T [v{\) +K<(f )T iv{_ 



(7) 

Vi (X) = 
For exam- 



ple, for the system with the potential ©, at the order u>°, 
where only the average of the potential over time is taken 
into account, V e ff = 0. The effective potential V e ff of 
is plotted in the inset of Fig. ^ The leading order 
(lj~ 2 ) is V^(x) = ^ 2 e~ 2 P x2 , which is clearly a dou- 
ble barrier. It is obvious that it traps the particle. This 
potential is always positive since it is the mean kinetic 
energy of the rapid oscillation around the slow motion, 
an energy which is coordinate dependent. For this it is 
instructive to note that V 2 = m£ 2 /2 in the leading order. 
In the order w -4 , on the other hand, the Hamiltonian 
cannot be expressed anymore in terms of an effective po- 
tential and terms that mix coordinates and momentum, 
like g(X)P 2 appear. Such terms result in corrections 
increasing with energy. It is also more complicated to 
understand terms of this order intuitively. 

Consider a quantum system with a Hamiltonian that is 
periodic in time, H(t + T) = H(t). Such a system can be 
treated using Floquet theory [ijj 0> EH • The symmetry 
with respect to discrete time translations implies that the 
solutions of the Schrodinger equation 
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arc linear combinations of functions of the form 



" U\{x,U)t) 



(8) 



(9) 



where A are the quasienergies and the corresponding 
quasienergy (or Floquet) states are u\[x, u>(t + T)) = 
u\(x,Ldt) with to — 2ir/T. This is the content of the 
Bloch-Floquet theorem in time. The states u\ are the 
eigenstates of the Floquet Hamiltonian 
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H. 



(10) 



These states have a natural separation into a "slow" 
part e~ % ~ (with the natural choice < X/h < to), which 
includes the information about the quasienergies, and a 
fast part u\(x,u>t) that depends only on the "fast" time 
t = tot. In the following an equation of motion for the 
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slow part of the dynamics is found as was done for clas- 
sical systems. It establishes a natural link between the 
separation into fast and slow motion in classical mechan- 
ics, that can be formalised by the theory of separation of 
time scales, and Bloch-Floquet theory in quantum me- 
chanics. For this purpose, following |l3j, we will look for 
a unitary gauge transformation e lF ^\ where F(t) is a 
Hermitian operator (function of x and p) defined at a cer- 
tain time t, which is a periodic function of time with the 
same period as H, such that in the new gauge the Hamil- 
tonian in the Schrodinger equation is time independent. 
In terms of the functions in the new gauge <f> = e ip, the 
Schrodinger Eq. JSJ is 



iTi 



d 



G4> 



where the Hamiltonian is 



G = e lF He~ lF 



ih 



de l 



-iF 



(ii) 



(12) 



Assume that such an operator F exists, so that G is 
time independent. Its eigenfunctions v\{x) evolve as 

4> x {t,x) = e- i %v x {x). (13) 

These states, in the original gauge, correspond to 

Va(*, x) = e- iF <f> x = e-^e- iF v x (x). (14) 

The function e~ lF v\ is periodic in time with the period 
of H and therefore ip\ of © is a Floquet state with 
quasienergy A (mod hui) where u\ — e~ lF v\. 

At high frequencies, F is found to be small, of the order 
of l/u. We expand G and F in powers of l/uo and choose 
F, so that G is time independent in any given order. The 
expansions are G — Xmlo G n /u n and F — XmLi Fn/u n . 
The calculation is performed by computing Gi in terms 
of Ft, ■ ■ ■ , Fi + i and then choosing F} + i so that Gi is time 
independent. The terms in (|12fl are calculated with the 
help of the operator expansion, 



°J F Be~ lF 



B + i 



F, B 



1 
2! 



F, 



F, B 



(15) 



where for the first term in (|12fl one takes B = H while 
for the second B = 

or 

The resulting time independent effective Hamiltonian 



V 2 - 1 



+ g(x)p 2 )+—V q + 0(u;- 5 ), 



(16) 



where V e ff(x) and g(x) are the classical terms (see (JSJ) 
and 0), while V q = ^(/ ( 



.(2), 



V, 



(3) 



is a quantum 



correction to the classical Hamiltonian that appears first 
in this order. The Hamiltonian (|16|) is presented in a 
form which is manifestly Hermitian. 

The effective Hamiltonian Ijl6|l is the main result of this 
work. Its classical limit is the classical effective Hamil- 
tonian © that can be obtained from JQ) by the canoni- 
cal transformation that is the classical limit of —TiF. It 
should be emphasized that in the derivation of (|16fl no 
semiclassical approximation was made. 

The perturbation theory that was developed here en- 
ables one to calculate not only the quasienergies that 
are the eigenvalues of G but also the corresponding 
quasienergy states. If the eigenfunctions of G are known, 
then the quasienergy (or Floquet) states can be computed 
up to order using equation (|14f) with 
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[Vi] 



(2)r 
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where the explicit expressions of terms of order u> 3 and 
cj~ 4 will be given elsewhere Q. 

The theory is demonstrated for the oscillating Gaus- 
sian J3Jl. The system demonstrates trapping by an os- 
cillating field, a phenomenon that is of physical interest. 
For this problem the effective potential 0, depicted in 
the inset of Fig. is a double barrier therefore it exhibits 
resonances. Each resonance is characterized by a com- 
plex energy E — iY /2. (For a relevant review see [lij.) 
For any resonance of (|16|l it is natural to look for the 
corresponding resonance of the original time dependent 
Hamiltonian 10. More precisely, one looks for the res- 
onances of the Floquet Hamiltonian (|10|) with H of 10. 
This is done numerically usin g a c ombination of the (t, t') 
method and complex scaling |26j . 

The energy Eq and the width To of the lowest (smallest 
real part Eg) quasienergy-resonance of J2J) are compared 
with the lowest resonance of the corresponding effective 
Hamiltonian (|16|) in Fig. ^ It is clear that for large fre- 
quencies there is excellent agreement. The results for the 
effective Hamiltonian truncated at orders lo~ a and uj~ 2 
are comparable at large frequencies. At low frequencies 
the results of order uu~ 2 turn out to be more accurate, in- 
dicating that the series do not converge at those frequen- 
cies. For comparison the characteristic frequency for the 
slow motion of the particle is in the range 0.1 — 0.3. 

It was found that the perturbation theory leads to a 
time independent effective Hamiltonian. This effective 
Hamiltonian may give physical insight, based on expe- 
rience with time independent systems, which is absent 
when examining the corresponding time dependent prob- 
lem. For example, consider a system where the time av- 
eraged potential consists of two barriers in addition to 
some high frequency time dependent perturbation (with 
vanishing average). If the perturbation is mainly in the 
region of the barriers one expects that the perturbation 
slightly raises the barriers in the effective Hamiltonian. 
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chaotic regions in phase space. This may hint that the 
perturbation theory may describe correctly only part of 
the phase space (whose fraction grows with ui). 

It is our great pleasure to thank M. V. Berry, N. David- 
son, N. Moiseyev, and V. Rom-Kedar for stimulating and 
inspiring discussions. We thank N. Moiseyev also for 
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This research was supported in part by the US-Israel Bi- 
national Science Foundation (BSF) and by the Minerva 
Center of Nonlinear Physics of Complex Systems. 



FIG. 1: The lowest quasienergy resonance of the oscillat- 
ing Gaussian 10, Eq (solid line) and To (dashed line), as 
a function of the driving frequency, compared to the lowest 
resonance of the effective Hamiltonian 1161 . Eq (diamonds) 
and To (circles), for 7 = 9 and /3 = 0.02 (in "atomic units" 
h — m — e — 1). Full symbols correspond to the effective 
Hamiltonian truncated at order uj~ 4 while empty symbols to 
order u> -2 . The effective potential (J7J is depicted in the inset, 
for lu = 1.5. 



In contrast, if the perturbation is in the region between 
the barriers it tends to raise the energy of the resonance. 
Therefore one expects that applying a time dependent 
perturbation in the region of the barriers will tend to 
increase the lifetime of the resonance {h/T) while apply- 
ing it in between the barriers will tend to decrease it. 
All the well developed techniques for time independent 
quantum systems can be used to compute the eigenvalues 
of G, in particular in the case where the eigenvalues and 

-2 

eigenstates of Go = frrr + Vb (x) are known. The effective 
Hamiltonian can also be used to predict trapping by os- 
cillating potentials that were so far investigated mainly 
numerically |27| . 

We have investigated the dynamics of high frequency 
driven general classical and quantum systems. High fre- 
quency perturbation theory, which exploits the idea of 
separation of time scales, was used to obtain an effective 
time independent Hamiltonian for the slow part of the 
classical and quantum motion. The spectrum of the effec- 
tive quantum Hamiltonian is the quasienergy spectrum 
of the time dependent system. This effective Hamilto- 
nian is computed to order w -4 in a perturbation theory. 
While in the order oj~ 2 the effect of the rapid oscillations 
around the slow motion could be expressed in terms of a 
classical scalar potential, the order w -4 involves both co- 
ordinates and momentum. Quantum corrections to the 
Hamiltonian also appear at this order. 

Some properties of the perturbation theory for F and 
G, such as its convergence and validity, are not well un- 
derstood and should be further studied. For instance, 
in one dimension the classical slow motion is integrable 
while the exact time dependent dynamics may have 
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